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ABSTRACT 

Extremely strong ionized Fe emission lines, with the equivalent width reaching about 

4000 eV, were discovered with ASCA from a few Galactic compact objects, including 

AX J2315-0592, RX J1802. 1+1804 and AX J1842.8-0423. These objects are thought 

^ . to be binary systems containing magnetized white dwarfs (WDs). A possible interpre- 

QQ ' tation of the strong Fe-K line is the line-photon coUimation in the WD accretion 

QQ I column, due to resonance scattering of line photons. The coUimation occurs when the 

^^ ■ accretion column has a flat shape, and the effect is augmented by the vertical ve- 

C"~~- ' locity gradient there, which reduces the resonant trapping of resonant photons along 

^D , the magnetic field lines. This effect was quantitatively confirmed with Monte-Carlo 

simulations. Furthermore, with ASCA observations of the polar V834 Centauri, this 

coUimation effect was clearly detected as a rotational modulation of the equivalent 

C^ . width of the Fe-K emission line. Extremely strong emission lines mentioned above 

^ i' can be explained by our interpretation consistently. Combing this effect with other 

' I X-ray information, the geometry and plasma parameters in the accretion column were 

;, . determined. 

^ ■ 
C/3 . Key words: radiation mechanisms: thermal - scattering - X-rays: stars ~ stars: white 

^ ' dwarfs - stars: individual: (V834 Cen) - methods: observation 

> 

Cu 1 INTRODUCTION lines, the implied abundances are usually sub-solar, such as 

OAtol solar for AM Her (Ishida et al. 1997), 0.63 ± 0.08 
A polar (or AM Her type object) is a binary system consist- gdar for EX Hya (Fujimoto, Ishida 1997), and ~ 0.4 solar 
ing of a low- mass main sequence star filling its Roche lobe foj. gg Cyg (Done, Osborne 1997). Therefore, we speculate 
and a magnetized white dwarf (WD) with lO'^ * G magnetic that the unusually high iron abundances of the present two 
field, which is strong enough to lock the WD spin with the wD binaries result from some mechanisms to enhance the 
orbital motion. Matter spilling over the Roche lobe of the li^e EW, rather than from high metallicities of the mass- 
companion star is captured by the magnetic field of the WD donating stars 

and accretes onto its magnetic poles, emitting hard X-rays Wg find a common feature in these two object, which 
via optically-thin thermal bremsstrahlung. may provide a clue to the strong iron K line. AX J2315-0592 
Extremely strong ionized iron emission lines have been exhibits a large (87 ± 2%, 57 ± 2%) orbital modulation in 
discovered from a few polars with the X-ray observation the 0.7 - 2.3 keV and 2.3 - 6.0 keV light curves, but al- 
by ASCA. For example, AX J2315— 0592 has a strong ion- most no modulation in the 6.0 - 10.0 keV band. Similarly, 
ized Fe-Kc Une centered at 6.84lQ;JgkeV, whose equivalent RX J1802. 1-1-1804 exhibits a large (~ 100%) orbital mod- 
width (EW) reaches 900t2oo ^^ (Misaki et. al 1996), and ulation amplitude below 0.5 keV with ROSAT observation 
RX J1802. 1-1-1804 has a strong Fe-K,, line with EW ~ 4000 (Greiner et al. 1998), but the ASCA light curves are ex- 
eV (Ishida et. al 1998). To interpret these strong line emis- tremely flat (Ishida et al. 1998). The lack of hard-band 
sions as thermal plasma emission, the plasma metallicity modulation implies that a constant fractional volume of the 
needs to much exceed one solar abundance; ~ 2 solar for accretion columns (which are optically thin to continuum 
AX J2315-0592, and 3.04±1. 47 solar for RXJ1802. 1-1-1804. X-rays) is observed throughout the rotational phase, and 
Although WD binaries often exhibit highly ionized Fe-K hence the inclination of the orbital plane i is rather small. 
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Figure 1. Schematic view of an accretion column on the WD in 
a polar system 



Then, the soft-band modulation must be due to changes in 
absorption by the pre-shock accretion flow, indicating that 
we observe down onto a single pole at the absorption max- 
imum. This in turn requires the magnetic co-latitude /? to 
be close to i. In short, these systems are inferred to have 
/3 ~ z ~ 0. In this paper, we call such a polars as "POLE 
(Pole-On Line Emitter)". 

We have another example of extremely strong iron K 
line emitter, an X-ray transient source AX J1842— 0423, 
discovered with ASCA in 1996 October on the Galactic 
plane in the Scutum arm region (Terada et al. 1999; here- 
after paper I). The most outstanding feature of this object 
revealed by the ASCA GIS is the very conspicuous emis- 
sion line at 6.781qj3 keV, whose EW is extremely large at 
4000^500° eV. To explain this line EW, a plasma metaUic- 
ity as high as 3.01'lQg solar abundance would be required. 
We found no periodicity over the period range of 62.5 ms 
to a few hours, the latter being a typical orbital period of 
polars. In view of the thin-thermal spectrum and the al- 
lowed source size of lo*-^"^^-^ cm, we have concluded in pa- 
per I that AX J1842-0423 is likely to another POLE, like 
AX J2315-0592 or RX J1802. 1+1804. 

The face-value metallicities of the three objects are so 
high that we regard these values to be unrealistic. Instead, 
we consider that the iron K line EW is much enhanced by 
some mechanism, which may be common to the POLEs. 
To account for the strong iron K line of three POLEs, the 
mechanism must account for line enhancement by a factor of 
3 or more. In this paper, we develop the possible explanation 
invoking resonance scattering (g3), which has been proposed 
briefly in paper I. In the present paper, we carry out Monte- 
Carlo simulations (93) to confirm the proposed mechanism, 
and verify the effect through ASCA observations of the polar 
V834 Cen (§|). 



2 LINE ENHANCEMENT DUE TO 
RESONANCE PHOTON BEAMING 

2.1 Geometrical beaming of Fe resonance line 

In a polar system, a flow of matter accreting onto each mag- 
netic pole of the WD is highly supersonic, so that a stand- 



ing shock is formed close to the WD. The matter is shock 
heated up to a temperature fc of a few tens keV [Appendix 
H equation (Al)]. The heated plasma cools by radiating 
Bremsstrahlung hard X-ray continuum and line photons, as 
it flows down the column, to form a hot accretion column of 
height h and radius r as illustrated in figure m We can as- 
sume that the ion temperature is equal to the electron tem- 
perature, because the ion to electron energy transfer time 
scale icq is much shorter tha n the coo ling time scale ttt.cooi 
[Appendix |A|, equations (A7) and (A6)]. 



The typical accretion rate of polars is M ~ lO^'' g s~^, r 
is typically 5 x lO'^ cm, and the velocity immediately beneath 
the shock front is typically u"^^ ~ 10* cm s"'^ [Appendix M 
equation (|A2|)], so that the electron density at the top of the 
hot accretion column is typically n'J' ~ 10^^ ^^ cm~^ [Ap- 
pendix kl equation (|A3|)]. The optical depth of the column 
for Thomson scattering is then given by 



Tt = 0.24 



V5 X lO'^cnJ 



7.7 X 1015cm" 



(1) 



and the optical depth of free-free absorption is much smaller 
[Appendix M equation ( A6 )] . Therefore, the column is opti- 
cally thin for both electron scattering and free-free absorp- 
tion, so the continuum X-rays are emitted isotropically. 

On the other hand, the optical de pth fo r the resonance 
scattering is calculated from equation ( AlC ) in Appendix M 
as 

^"^ ^ ^^ V7.7X 1015cm- V U.68X 10-^ J U x lO^cmJ ' ^^^ 

at the energy of the hydrogenic iron K^ line, where Ap„ is 
the abundance of iron by number, which is normalized to 
the value of one solar. Thus, the accretion column is opti- 
cally thick for resonance lines, and the resonance line pho- 
tons can only escape from positions close to the surface of 
accretion column. If the accretion column has a flat coin- 
shaped geometry, and our line-of-sight is nearly pole-on to 
it, we will observe the enhanced Fe-K lines. We call this ef- 
fect "geometrical beaming". However, this effect can explain 
the Fe-K line enhancement up to a factor of 2.0 [Appendix 
^, which is insufficient to explain the factor ~ 3 enhance- 
ments observed from POLEs. An additional mechanism is 
clearly needed. 

2.2 Additional collimation by vertical velocity 
gradient 

In the accretion column of a polar, both temperature kT 
and u decrease, and rie increases from the shock front to- 
ward the WD surface. Numerically, the vertical profiles of 
these quantities are calculated by Aizu (1973), as a func- 
tion of vertical distance z from the WD surface, as equation 
( A4 ) in Appendix |A[ Because of Doppler shift caused by this 
strong vertical velocity gradient in the post-shock flow, the 
resonance line energy changes continuously in the vertical 
direction. Let us consider, for example, that a line photon 
of rest-frame energy Eg is produced near the bottom of ac- 
cretion column, z ~ O.lh, where the emissivity of He-like 
iron line photon becomes maximum, and that this photon 
moves vertically by its mean free path of resonance scatter- 
ing Zr which is given as 



Zr 



(anAFaTle 
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photon produced at zlh = 0.1 




Figure 2. (left) The cumulative shift of the resonance energy (solid line) for upward-moving Fe-K line photons, which are produced at 
2 = O.lh, 0.2h, 0.3h, QAh, and Q.5h. The thermal width of the resonant Fe-K line is given by the dashed curve. Separations of the adjacent 
dotted vertical lines specify the mean free path Zr of resonant photons, (right) Cross section of resonance scattering for an iron Kq line 
photon, which is produced at z/h = 0.10 and moves vertically. The resonance centroid energy corresponds to AE^^ shown in the left 
panel, and the resonance energy width corresponds to AS*''. 



= 1.3 X 10' 



U.68 X 10-5 J 



7.7 X 1015cm 



(3) 



Because h is given as equation (A5) [Appendix |A|, l^ is 
roughly equal to 0.07/i(|-)^'^. Over this distance along the 
z-direction, u will chang e by Am = 4.0 x 10^ cm s~^ at 
z ~ 0.1 /i from equation (A4), so the resonance energy for 
the line photon shifts due to Doppler effect over the same 
distance by 



AS" 



£^0 



Au 



0.9 



V6.695keVy 



0.9 X 10«cm 



eV, 



(4) 



where c is the light velocity. Figure shows the change 
of /S.E'°^ for upward-moving photons produced at various 
depths of the accretion column. 

The width A_ E of the resonance scattering [see Ap- 
pendix |A| equation ( A9)] is determined by the natural width 
of line photon (~ 1 eV for iron K ion) and by the thermal 
Doppler broadening AiJ"*. Because the thermal velocity of 
ion of mass rrii reaches 



kT""- 



2.4 X 10 



'©* 



C. 



TTli 



-h ( kr'^ \ ' 



V56mH/ ylGkeVy 



cm s ^{5) 



with run the mass of a hydrogen atom, the thermal Doppler 
width of the resonance becomes 



AE'^iz) 



Eo 



t'i 



V6.695keV/ \l6keV I 



V 56mH / 



(6) 



Here we normalized the equation to the iron atom, mi = 
56mH. The dashed line in figure H shows the position depen- 



dence of AiJ"', below which the resonance scattering does 
occur. 

Figure H clearly shows that, if a line photon grad- 
ually moves upward through repeated scattering, its en- 
ergy becomes different from the local resonance energy, by 
an amount AE^'^ [equation (H)] which eventually becomes 
larger than the thermal width Ai?"" [equation (0)]. Numer- 
ically, this ratio for a photon produced ai z — zo can be 
described as 



AE^'-jz) 
AE^^{z) 



Am 



«=''|(2//l)*-(zo//l)*} 



= 5.5 



2.2 



V56mH/ V 0.615/ 
V56mH/ V 0.615/ 



(7) 



(8) 



where 5z = z — zq, and ^ is the mean molecular weight 
(/i = 0.615 for a plasma of one-solar abundance), and we 
used the relation of kT^^ = S^mn • {u"^)'^ from equations 
( |A1| ) and ([A2|). Then, the photon is no longer scattered ef- 
ficiently, and can escape out. This effect does not occur in 
the horizontal direction because of little velocity gradient. 
As a result, a resonant line photon produced near bottom of 
the accretion column thus escapes with a higher probability 
when its net displacement due to random walk is directed 
upward, rather than horizontal. In other words, iron K line 
photons are coUimated to the vertical direction. We call this 
effect "velocity gradient beaming" . The essence of this effect 
is that the mass of iron, mi in equation nm, is heavy enough 
for the bulk velocity gradient to overcome the thermal line 
broadening. 
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Figure 4. Temperature dependence of the Fe line emissivity 
(Mewe et al. 1985) adopted in the calculation. 



Figure 3. Definition of the cylinder in our simulation. 



3 MONTE-CARLO SIMULATION 

The collimation effect described in section y involves signif- 
icant scattering processes, which make analytic calculation 
difficult. In this section, we accordingly examine the pro- 
posed effect using Monte-Carlo simulations. 



3.1 overview of the simulation 

Consider a simple cylinder of height h and radius r, filled 
with X-ray emitting plasma. We describe the vertical de- 
pendence of fcr, Me, and u in the cylinder by equation ( A4). 
We then produce a number of Monte-Carlo iron line photons 
in proportion to its emissivity, which in turn is determined 
by kT and n^ there. We isotropically randomize their ini- 
tial direction of propagation in the rest frame of iron nuclei. 
The line energy of each iron photon is also randomized; i.e. 
the average value of the energy is Doppler shifted in the 
frame of observers according to the bulk velocity law [eq. 
(y)], and its dispersion is determined by the local thermal 
motion [eq. (H)]. We trace the propagation of each line pho- 
ton with a constant step length, which is taken to be 1/100 
of the mean free path of resonance scattering at the bottom 
of the cylinder, where the temperature falls below 1 keV. At 
each step, the behavior of photon is determined by the cal- 
culated probabilities of resonance scattering and Compton 
scattering. We follow the propagation of each photon until 
it moves outside the cylinder. 



3.2 basic condition 

The temperature dependence of iron line emissivity in op- 
tically thin plasma has been calculated by many authors, 
and we here adopt the calculation by Mewe et al. (1985) as 
shown in figure \ We consider four species of iron K line 
photons; those of H-like resonance Kc line (6.965 keV), He- 
like resonance K^ line (6.698 keV), He-like intercombination 
line (6.673 keV), and He-hke forbidden line (6.634 keV). The 
emissivity of iron line photons per unit volume, in erg s~^ 
cm~^, is described for each species as 



where P'{kT) in erg s~^ cm~^ is the value shown in fig- 
ure M. The iron density npo, in cm"'', is calculated assuming 
one solar abundance. The position dependence of iron line 
emissivity is determined, through e qua tion (H), by the z de- 
pendence of kT and npo [equation (A4)]. 

Treatment of the resonance scattering process, taking 
into account both the bulk flow and the thermal motion of 
ions, is a key point of the present simulation. For each line 
photon being traced, its scattering probability at the i-th 
step fi is calculated as oc npc(fi)crRs(-Ei"), where nFo(ri) is 
the local Fe-ion density at the position ri [see Appendix H 
equations (A3) and (A4)], ctrs is the cross section for the 
resonance scattering given by eq. (|A9|) [Appendix H, and 
i?^" is the Doppler-shifted energy of the incoming photon 
measured in the rest frame of a representative Fe-ion at ri. 
The velocity Wi of this Fe-ion, relative to the observer, is 
expressed as a sum of the bulk flow velocity m at r^, and a 
random thermal velocity v. We specifically calculate as 



{1 + {wi-i — Wi) ■ ei"i^i/c} , 



(10) 



P = P'{kT) 



TlcTlpci 



(9) 



where Wi-i is the observer-frame velocity (bulk plus ran- 
dom) of the Fe-ion that scattered the line photon last time, 
E°^\ is the outgoing photon energy as expressed in the rest 
frame of the previous scatterer, eiJ'i^i is the unit vector along 
the photon propagation direction from the (i — l)-th to the 
i-th scattering site [figure H. If the scattering occurs at ft, 
we randomize the line photon energy from i5|" to £°"' ac- 
cording to the natural width, and isotropically randomize 
the direction of the outgoing photon, both in the rest frame 
of the present scatterer. If, instead, the scattering does not 
occur at fi, we preceed to the next step witouht changing its 
direction. Thus, our calculation autonatically incUudes both 
the bulk-flow Doppler effect and the thermal broadening. 
However, we do not consider energy shifts by the ion recoil, 
which is completely negligble. The scattering probability for 
a non-resonant photon is set to 0. 

In addition to the resonance scattering, we must con- 
sider the Compton scattering process; the energy E is shifted 
to E/{l-\ — ^ (1 — cos6'cMp)}, where ^cmp is the Compton 
scattering angle. For an iron Kq photon with energy E ~ 6.8 
keV, a Compton scattering with 6cmp > 10° will change the 
photon energy beyond the resonance energy width of a few 
eV: then the resonance scattering can no longer occur after 
a large-angle Compton scattering. We take this effect into 
account in our simulation, using the probability distribu- 
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Figure 5. The volume emissivity of Fe-K line photons in the assumed accretion column, shown as a function of zlh. The plasma 
temperature is kT^^ = 16 keV (left panel), and 60 keV(right panel). 
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Figure 6. Schematic illustration of each step in our calculation; 
fi is the calculation site of i-th step, Wi is the ion motion there 
in the observer's frame, which is a sum of bulk flow motion and 
random thermal motion, and Ei is the photon energy in the rest 
frame of i-th ion (See the text). 



tion of ^CMP by the Klein-Nishina formula, which is almost 
identical to the classical formula for the energy of iron lines. 
The differential scattering cross section and the energies of 
scattered photons are calculated in the rest frame of the 
currently scattering electron, so that the anisotropic effects 
caused by the bulk motion of electron is also included. Note 
that we neglect the process that the energy of a Compton- 
scattered continuum photon comes accidentally into the res- 
onance energy range, since we do not generate continuum 
photons in the Monta-Carlo simulation. 

3.3 results 

First we simulated the simplest case wherein the plasma is 
hydrostatic with a single temperature and a single density: 
i.e. u is set to and there is no vertical gradient in kT or 
rie. The angular distributions of line photon flux, calculated 
under this simple condition for various densities, are shown 
in figure |7[ The results confirm that the photons are emit- 
ted isotropically when the plasma density is low, but as the 
density increases, the geometrical beaming becomes progres- 




10 20 30 40 50 60 70 an 90 

polar aQgle (degrees) 

Figure 7. Angular distributions of resonant line photons emer- 
gent from a simple thin thermal plasma, simulated by neglecting 
the vertical motion (u = 0) and neglecting the vertical gradient 
in n and kT. Abscissa is the angle 9 defined in figure pi and ordi- 
nate is the photon flux per unit steradian normalized to the value 
of an isotropic emission. The plasma parameters are set to the 
nominal values; kT = 16 keV, r = 5 X 10^ cm, and h = 1.9 X 10^ 
cm. The solid line represents the analytic solution when only the 
column surface shines (equation Blj) . The crosses show the calcu- 
lated results for various densities as given in the figure. 



result agrees nicely with the analytic solution which assumes 
a completely optically-thick condition; i.e. line photons emit 
only from the surface of the accretion column [equation (Bl); 
see Appendix H. This verifies proper performance of our 
Monte-Carlo simulation. 



Next we have fully considered the vertical gradient in 
kT,u, and n^ [eq. (A4)]. Figure H shows the calculated an- 
gular distribution of He-like iron line when the relevant pa- 
rameters are set to the nominal values in equations (|Al|), 
(|A2|) , and ( A.3 ) . The resonance line fiux is thus enhanced in 



sively prominent. At We 



10^ 



cm , the Monte-Carlo 



the vertical direction more strongly than in figure m This 
reconfirms the physical beaming mechanism we proposed 
in § y. The Compton scattering is confirmed to reduce the 
collimation only slightly. Finally, the intercombination pho- 
tons, which are free from the resonance scattering, exhibit 
a nearly isotropic distribution. 
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Figure 8. The same as figure M, but calculated for the He-like 
iron Kq lines considering the vertical structure of the accretion 
column. The plasma parameters are set to the nominal values 
given in appendix ^: kT''^ = 16 keV, u''^ = 0.9 X 10* cm s"^, 
njii = 7.7 X 10^5 cm-3, r = 5 x lO'' cm, and h = 1.9 x 10^ cm. 
Thick solid crosses represent the results considering the Comp- 
ton scattering, while dashed crosses represent those neglecting 
the Compton process. The dotted crosses show the profile of the 
intercombination line. The thin crosses show the same profile as 
presented in figure M (nominal case). 



We repeated the Monte-Carlo simulations by changing 
h/r, n"^, kT, and u"^, around their baseline values of kT — 
16.0 keV, u"'' = 0.9 x 10* cm s-\ nf = 7.7 x 10^^ cm-^ 
r = 5 x 10^ cm, and h = 1.9 x 10^ cm (Appendix 0). Figure 
n(j summarizes the obtained results in terms of the beaming 
factor 



c= 



/(o) 



/o^ md 



cos 9 



(11) 



where f{d) is the angular distribution of line photon flux, 
such as is shown in figure B The beaming factor ^„ increases 
as h/r decreases (i.e. coin shaped column), or density in- 
creases (figures |l0p , and b) . However, when the density ex- 
ceeds ~ 10^® cm ~^, the beaming effect diminishes again, 
because of large-angle Compton scattering. This inference 
is achieved by comparing results with and without Comp- 
ton process (figure Mb). The WD mass dependence of 5„, is 
small (figure 0): it increases slightly with mass increases be- 
cause shock velocity increases with deeper gravity potential, 
and it starts decreasing because density decreases. These re- 
sults show clearly that the strong coUimation of He-like iron 
Kq photons, with ^„ > 2, is possible under reasonable con- 
ditions. 



4 AN OBSERVATIONAL APPROACH 

In order to experimentally verify our interpretation, it is 
necessary to measure the equivalent width of the resonant 
Fe-K lines as a function of viewing angle. For that purpose, 
we may utilize a polar of which our line-of-sight relative to 
the magnetic axis changes from ~ 0° (pole-on) to ~ 90° 
(side-on), as the WD rotates. Among the polars with well 
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Figure 9. The beaming factor 5„ with various WD mass. The 
H-like iron line is represented by solid lines and the He-like iron 
resonance line by dashed lines. The baseline condition is given in 
the text, and the temperatures and the densiti es c hange wit h WD 
mass and WE radius as shown in equations ( |Al| ) and ( |A2| ). The 
accretion rate to WD M for calculation is set to 1 X ICr^ g s~^ 
and 5 X 10^^ gs"!. 



determined system geometry (via optical, ultraviolet and in- 
frared observations), V834 Cen is particularly suited: its or- 
bital plane is inclined to our line-of-sight by i = 45±9°, and 
its magnetic co-latitude is /3 = 25 ± 5° (Cropper 1990). As 
a result, our line-of-sight to the accretion column changes 
from 20° to 70° . 



4.1 observations of V834 Cen with ASCA 

ASCA has four X-ray Telescopes (XRT:Selemitsos et.al 
1995), and its common focal plane is equipped with two 
Gas Imaging Spectrometers (GIS: Ohashi et al. 1996; Mak- 
ishima et al. 1996) and two Solid-state Imaging Spectrom- 
eters (SIS: Burke et al. 1991; Yamashita et al. 1997). The 
ASCA observation of V834 Cen was carried out for about 20 
ksec from 1994 March 3.63 to 1994 March 4.13 (UT), and 
about 60 ksec from 1999 February 9.93 to 1999 February 
10.72 (UT). In these observations, the GIS was operated in 
PH-nominal mode, which yields 0.7-10.0 keV X-ray spec- 
tra in 1024 channels, and the SIS was operated in 1-CCD 
FAINT mode, which produces 0.4-10.0 keV spectra in 4096 
channels. The target was detected with a mean count rate 
of 0.171 c s"^ per GIS detector and 0.259 c s~^ par SIS de- 
tector in 1994. The corresponding count rates were 0.191 c 
s"^ and 0.253 c s"^ in 1999. 

For extracting the source photons, we accumulated the 
GIS and SIS events within a circle of radius 4'5 centered on 
V834 Cen, employing the following data-selection criteria. 
We discarded the data during the ASCA passing through the 
South Atlantic Anomaly. We rejected the events acquired 
when the field of view (FOV) of ASCA was within 5° of the 
Earth's rim. Furthermore for the SIS, we discarded the data 
acquired when the FOV is within 10° of the bright Earth 
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Figure 10. A summary of Monte-Carlo simulations. The beaming factor ^„, is presented for various column shape and density. The 
H-like iron line is represented by solid lines and the He-like iron resonance line by dashed lines. The baseline condition of the calculation 
is given in the text, a) The shape dependence of §„, where the volume r^h is fixed. Calculation is performed for three densities as 
specified in the figure, b) The density dependence of ^m for three different column shapes (coin shape, nominal, cylinder shape). The 
coin shaped column has r = 1 X 10* cm and h = 1 X lO'^ cm, while the cylinder shaped column has r = 1 X 10^ cm and h = 1 X 10* cm. 
The beaming factors without Compton scattering are also shown. 



rim and those acquired near the day-night-transition of the 
spacecraft. 



4.2 Light curves 

Figure O shows the energy resolved Ught curve of V834 
Cen obtained with ASCA SIS-I-GIS folded by its rotational 
period, 1.69194 hr (Schwope et al. 1993). The phase is co- 
herent between the two light curves. The pole-on phase of 
V834 Cen is determined by the optical photometry and po- 
larimetry as shown in table 1 of Bailey et al. (1983), which 
corresponds to a phase (j) ~ 0.65-0.85 in our X-ray light 
curve. We can recognize small dips in light curves of the 
softer two energy bands in figure [111 at ~ 0.65 in 1994 
and ~ 0.86 in 1999. These dips are thought to arise from 
photoelectric absorption by the pre-shock matter on the ac- 
cretion column. Therefore, the pole-on phase is consistent 
between the optical and X-ray datasets. The folded light 
curve in the iron line energy band (6.2 - 7.2 keV) exhibits 
a hump at or near this pole-on phase, suggesting that the 
proposed line photon enhancement is indeed taking place in 
this system. However, detailed examination of the modula- 
tion of iron line emission needs phase-resolved spectroscopy, 
performed in the next subsection. 



4.3 spectral analysis 

Because the folded light curves (figure hll) have different 
shapes between the two observations, and because spectral 
information of GIS-3 in 1994 was degraded by a tempo- 
rary malfunctioning in the GIS onboard electronics, we use 
only the 1999 data for spectral analysis. We have accumu- 
lated the GIS (GIS2 + GIS3) and SIS (SISO + SISl) data 
over the pole-on phase (<j!> = 0.73 ± 0.25) and side-on phase 
{(p = 0.23 ±0.25) separately. We subtracted the background 
spectrum, prepared by using the blank sky data of the GIS 
and SIS. The spectra, thus derived and shown in figure Pl2l 
exhibit an absorbed continuum with strong iron Kq emission 
lines over the 6.0-7.2 keV bandpass. 



4-. 3.1 Continuum spectra 

We attempted to quantify the 0.8 - 10.0 keV continua, ne- 
glecting for the moment the line energy band of 6.0-7.2 keV. 
However, the simplest spectral model for polars, namely 
a single temperature bremsstrahlung continuum absorbed 
by one single column density, failed to reproduce either of 
the observed spectra (figure |l^ left; Model 1 in table ??). 
This failure is not surprising, considering that polars gen- 
erally exhibit multi-temperature hard X-ray emission with 
complex absorption by the pre-shock absorber (Norton and 
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Figure 11. Energy resolved light curves of V834 Cen obtained with ASCA, folded at 1.69194 hr. Phase = corresponds to HJD 
2445048.9500, which is common to the 1994 and 1999 light curves. Each light curves is shown for two cycles. Left panel shows the 
observation in 1994, right panel in 1999. 



Pole on spectrum 



Side on spectrum 




Energy (keV) Energy Q<eV) 

Figure 12. ASCA X-ray spectra of V834 Cen taken in 1999, shown without removing the instrumental responses. Left panel shows the 
spectra in the pole-on phase (</> = 0.73 ± 0.25), and right panel those in the side-on phase (</> = 0.23 ± 0.25). The CIS and SIS data are 
shown with crosses. The solid crosses show the best-fit model consisting of single temperature bremsstrahlung continuum with partially 
covered absorption and three narrow Gaussians (see text). The best-fit parameters are shown in table ?? Model 4 and in table B^ (full 
band fitting). Lower panels show the fit residuals. 



Watson 1989). The extremely high temperature obtained 
by this simple fitting, > 200 keV, is presumably an artifact, 
compared with the temperature of 14.7 keV measured with 
Ginga. 

To accurately estimate the hottest component of the 
continuum, avoiding complex absorption in soft energies, we 
then restricted the fit energy band to a narrower hard en- 



ergy band of 4.5 - 10.0 keV (Model 2 in table ??). This 
lower limit (4.5 keV) was determined in a way described 
in §3.2 by Ezuka and Ishida (1994). In this case, A^h is de- 
termined solely by the depth of the iron K-edge absorp- 
tion at ^ 7.1 keV. This model has been successful on the 
spectra of both phases, yielding a temperature consistent 
with the Ginga value. Considering that the K-edge absorp- 
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si ngle kT, single A^h 



m ulti fcT, single A^h 




Figure 13. The same spectra as in figure [13 left (pole-on). Left panel shows the best-fit model of single temperature bremsstrahlung 
absorbed by a single column density. The best-fit parameters are shown in table ?? Model 1. Right panel shows the best-fit model of 
double temperature bremsstrahlung absorbed by a single column density, corresponding to table ?? Model 3 (free kT). 



a. s ingle narrow 



b. single broad 



c. double narrovy 
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Figure 14. Phase averaged SIS spectra in the iron K-line energy band, fitted with a single narrow Gaussian (left), a single broad 
Gaussian (center), and a double narrow Gaussian model (right). The continuum spectrum is represented by a single temperature and a 
single column density in the 4.5 — 10.0 koV band (narrow band fitting). The best-fit parameters are shown in table Q. 



tion in the observed spectra are relatively shallow, the value 
of A^H obtained in this way is thought to approximate the 
covering-fraction-weighted mean value of multi-valued ab- 
sorption. The mean value of A^h in the accretion column is 
hence inferred to be ~ 10^'^ cm~^. 

We next fitted the original 0.8 - 10.0 keV spectra (but 
excluding the iron K line region) by a two-temperature 
bremsstrahlung with single A''h, and obtained acceptable re- 
sults (right panel of figure [13^ and Model 3 in table ??). How- 
ever, the first temperature kTi is still poorly determined; 
the fit became unacceptable for the pole-on spectra when 
we fixed kTi to the Ginga value. Furthermore, the obtained 
A'h is about 1 x 10^^ cm~^ in either case, which is not consis- 
tent with the inference from the narrow band fitting. Thus, 
we regard Model 3 as inappropriate. 

A fourth spectral model we employed consists of a sin- 
gle temperature bremsstrahlung and double-valued photo- 
electric absorption (A''hi and Nh2), which is so-called par- 
tial covered absorption model (PCA model) . This model has 
been fully acceptable for both phases (Model 4 in table ??). 
The obtained A^h2 is consistent with that suggested by the 
narrow band fitting, and the fit remained good even when 
we fix the temperature to the Ginga value. We therefore uti- 
lize this model (Model 4) as the best representation of the 
continua for both phases. The solid curves in figure p_2 refer 
to this modeling. 



4.3.2 The Fe-K lines 

Having quantified the continuum spectra of V834 Cen, we 
proceed to the study of the iron K-line. For this purpose, we 
employ the phase-averaged spectrum, and again limit the 
energy range to the 4.5 - 10.0 keV band to avoid complex ab- 
sorption structure in lower energies. As a first-cut attempt, 
we modeled the iron line with a single Gaussian model, 
while represented the continuum with a single-temperature 
bremsstrahlung absorbed with a single column density (to 
reproduce the iron K edge), but the fit failed to reproduce 
the line profile as shown in figure hja. A broad Gaussian 
model with g ~ 0.28 keV has been found to be successful 
(table Gt figure [l4|b), but the obtained line centroid energy 
is too low for ionized Fe-K species that are expected for a 
plasma of temperature ~ 10 keV. We can alternatively fit 
the data successfully with two narrow Gaussians (table Et 
figure |l4|c) , where the centroid energy of the first Gaussian 
turned out to be consistent with that of the fluorescent iron 
Kq fine (6.40 keV); that of the second Gaussian comes in 
between those of He-like iron Kq line (6.65 - 6.70 keV) and 
H-like line (6.97 keV). This means that the second Gaussian 
is in reality a composite of the H-like and He-like lines. 

We therefore employed a line model consisting of three 
narrow Gaussian, each having a free centroid energy and a 
free normalization. We have then obtained an acceptable fit, 
with the three centroid energies consistent with the Fe-K line 
energies of the neutral. He-like, and hydrogen-like spices, as 
shown in figure Il5l These three lines have been observed in 
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Table 1 . Best Fit continuum parameters for the GIS and SIS spectra of V834 Cen" 



Model'' 


fcTi 


fcTa 


A^Hl 


Nh2 


Gov. Frac= 


Xl (dof) 




keV 


keV 


Xl022 cm-2 


Xl022 cm-2 


% 




POLE_ON 


phase 












Model 1 


>200 


- 


n 19+002 
^•1^-0. 02 


- 


- 


2.19 (188) 


Model 2'' 


lO.Stfg* 


- 


14 5+8-1 


- 


- 


0.49 (94) 


Model 3 


> 158 


n 1Q+001 


n 09+O.O8 


- 


- 


1.01 (186) 




14.7= 


n ic:+0.01 


1 91 +0.08 
^•^^-0.08 


- 


- 


1.64 (187) 


Model 4 


> 164 


- 


< 0.02 


741+2.26 


41.31^;^ 


0.91 (186) 




14.7= 


- 


< 0.08 


QQ + l-SO 
»-»'3_i.59 


60.41^;^ 


1.00 (187) 


SIDE_ON phase 












Model 1 


> 171 


- 


< 0.04 


- 


- 


1.00 (174) 


Model 2^* 


12.4+1*0 


- 


13.9+^1-^ 


- 


- 


0.47 (70) 


Model 3 


>200 


PI r,c.+0.03 


0.36«:}^ 


- 


- 


0.85 (172) 




14.7= 


n 1O+0.01 


0.78«;1« 


- 


- 


1.135 (173) 


Model 4 


>49.3 


- 


< 0.02 


12.7+^^5 


35.6tlH 


0.81 (172) 




14.7= 


- 


< 0.06 


11 2+'l-5 


50.7tli 


0.85 (173) 


Phase average 












Model 2<* 


11 7+15.2 
ll-'-4.4 


- 


u.otli 


- 


- 


0.49 (122) 


Model 4 


14.7= 


- 


< 0.06 


9.79ll:«9 


56.9l2;i 


1.07 (275) 



Excluding the Fe Ka line band (6.0 - 7.2 keV). 

Model 1/2: single Nh, single kT. Model 3 : single A^h, multi kT. 

Model 4: multi A^h, single kT. 

The covering fraction (%) of A^hi. 

Fitting in the 4.5 — 10.0 keV band. For other models, the 0.8 - 10.0 keV band is used. 

Gontinumn temparature fixed at the value measured with Ginga (Ishida et al. 1991). 




6.4 6.6 6.8 7 

Line Center Energy [keV] 



7.4 



this model. For consistency check, we have expanded the fit 
energy range back to 0.7 - 10.0 keV, and performed full-band 
fitting, employing the continuum Model 4 and the three nar- 
row Gaussian for the iron K-line. The results, presented in 
table H, are generally consistent with those from the narrow- 
band analysis. 

The temperature and A^'h obtained in these final fits are 
thus the same between the two phases within errors, and 
the EW of fiuorescent and H-like lines are also consistent 
with being unmodulated. In contrast, the EW of the He-like 
line is enhanced by (^poic/Csido)"'"' = 1.871q'54 times (narrow 
band fitting) or 1.9llQg4 times (full band fitting) in the 
pole-on phase compared to the side-on phase, with > 90 % 
statistically significance. 



Figure 15. The best-fit line centroid energies with 68 % errors in 
terms of the three narrow Gaussian model. From the top to bot- 
tom panels, the results for the phase average, pole-on phase, and 
side-on phase are shown. The upper plot in each panel shows the 
results when three centroid energies are left free. The lower plot 
shows the results when the centroid energy of the first Gaussian 
is fixed at 6.40 keV. The fit energy range is 4.5 - 10.0 keV. 



the spectra of many polars with ASCA (Ezuka and Ishida 
1999). Hereafter, we adopt the three narrow Gaussian model 
in quantifying the iron lines of V834 Cen. 

As our final analysis, we have repeated the three- 
Gaussian fitting to the spectra, by fixing the line centroid 
energy of the first Gaussian (identified to a fiuorescent line) 
at 6.40 keV. The result is of course successful, and the ob- 
tained parameters are given in table |3| as well as figure |l5| 
Figure Iq shows the phase-resolved SIS spectra fitted with 



5 DISCUSSION 

5.1 Origin of the line intensity modulation in 
V834 Cen 

In order to verify our interpretation, mentioned in section 
H we observed the polar V834 Cen with ASCA. The EW of 
the He-like iron-Ka line has been confirmed to be enhanced 
by a factor of (^poie/Csido)°''° ~ 1-87 ± 0.54 in the pole-on 
phase relative to the side-on phase. Can we explain the ob- 
served rotational modulation of the He-like iron line EW 
of V834 Cen by some conventional mechanisms? An obvi- 
ous possibility is that the temperature variation modulates 
the iron line intensity. However, the continuum temperature 
determined by the narrow band fitting (table w) is almost 
constant within the error to explain the rotational modu- 
lation of the line intensity ratio. Alternatively, the bottom 



© 2001 RAS, MNRAS 000, Q-^ 



X-Ray Beaming due to Resonance Scattering in the Accretion Column of Magnetic Cataclysmic Variables 11 



Table 2. Best Fit parameters of the iron Kq line" 



Line model 



iron K„ line 



statistics 



L.C. l*" 
(keV) 



CTl 

(keV) 



L.C. 2'' 
(keV) 



0-2 

(keV) 



X5 (dof) 



single narrow 
single broad 
double narrow 



6.63 

6 

6.45 



+0.04 
0.05 
+0.05 
-0.05 
+0.08 
0.11 



0.0001= 

'^■'^='-0.05 

0.0001^= 



6.79 



+0.23 



0.0001= 



1.37 (69) 
0.66 (68) 
0.91 (67) 



The phase averaged spectrum. Only SIS data are used for the fitting. 

Fitted with a single temparature and single column density model in 4.5-10.0 keV. 

Fixed the continuum temperature and A'^h to the Model2 in table ?? 

Line center energies. 

Fixed. 



Pole-on spectrum 




iii*yw^^^^ 



Side-on spectrum 




}ttWVMfW¥V+ 



Figure 16. Phase resolved spectra over the Fe-K line energies, fitted with three narrow Gaussians. The central energy of the first 
Gaussian is fixed at 6.4keV. The best-fit parameters are shown in table H. 



part of the accretion column, where the emissivity of the 
He- like iron line is highest (figure m, may be eclipsed by 
the WD surface so that the He-like line is reduced in the 
side-on phase. However, such a condition does not occur un- 
der the geometrical parameters of V834 Cen. Furthermore, 
we do not see any dip in the rotation-folded hard X-ray 
light curve. Therefore, this explanation is not hkely, either. 
A third possibility is that a part of iron line photons come 
from some regions other than the accretion column. In fact. 
X-rays from polars are known to be contributed by pho- 
tons reflected from the WD surface (Beardmore et al. 1995; 
Done & Beardmore 1995; Done & Magdziarz 1998). This 
mechanism can explain the production of fluorescent lines, 
and possibly its rotational modulation, but cannot produce 
highly ionized iron lines. From these considerations, we con- 
clude that the observed iron line modulation is difficult to 
account for without appealing to the resonance scattering 
effects. 

Our next task is to examine whether the geometrical 
coUimation mechanism (§ 2.1) can explain the observation. 



When the optical depth of resonance scattering is very high 
and hence the line photons come solely from the surface of 
the accretion column, the angular distribution of the iron 
line intensity is given analytically by equation (Bl) in Ap- 
pendix |b|. By averaging this distribution under the geometry 
of V834 Cen, and taking into account the exposure for the 
two phases, we have calculated the expected line beaming 
factor ^ for the ASCA data as shown in Figure JVA Thus, 
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to explain the observed enhancement of He-like line solely 
by the geometrical beaming, a very fiat coin-shaped column 
with h/r ~ 0.2 would be required. Furthermore, the ob- 
served ^ should be lower by about 30-40% than the ideal 
calculation in Figure [ITI because finite optical depths, ex- 
pected under a typical plasma density of polar accretion 
column, reduce the enhancement (Figure Vm . We hence con- 
clude that the geometrical mechanism alone is insufficient to 
explain the observed He-like iron line enhancement in V834 
Cen (§ 0), and hence the additional collimation due to the 
velocity gradient effect is needed. 

Then, are all the observed results consistently ex- 
plained in our picture that incorporates the geometrical and 
velocity-gradient effects? At the temperature of V834 Cen 
(~ 15 keV), the resonant photons in fact contribute about 
65% to the observed He-like iron Ko, line, the rest coming 
from intercombination (~ 20%) and forbidden (~ 15%) lines 
which are free from the resonance scattering effects (figure 
kl). As a result, the line collimation is expected to be some- 
what weakened in the ASCA spectra, where we cannot sepa- 
rate these unmodulated lines from the resonance line. After 
correcting for this reduction, the true value of ^ for the He- 
like resonance line is calculated to be ~ 2.3. This is within 
the range that can be explained by our scenario. 

How about the H-like iron Kc line? Even though it con- 
sists entirely of resonant photons, its modulation in the V834 
Cen data has been insignificant (table H). Presumably, this 
is mainly due to technical difficulties in detecting this weak 
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Table 3. Best Fit parameters of the iron Kq line with three narrow gaussans. 



Phase" continuum 


Fluorescent 


He-like 


H-like 


X5 (dof) 


kT 


Nh 


EW 


I.e." 


EW 


I.e." 


EW 




(kcV) 


Xl022 cm-2 


(eV) 


(keV) 


(eV) 


(keV) 


(eV) 




Narrow Band Fitti 


ng'^ 














Average 9.3+^°9'' 
Pole-on 9.6J:fgf 
Side-on 9.0tfo° 


15.4+^8 

i5.6l«:° 

1fi R+13-4 


244+^3 
242l?^2 
270tl- 


6 70+°°" 

"•'"-0.04 
fi CQ+0.07 

c cn+0.70 


3521^; 
417t;°^ 
223t?^i 


7 07+0'08 
'•"'-0.44 

7 04+°-" 

7 14+0. 16 
'•l"t-0.48 


192+97 
990+123 

2231^^ 


0.67 (154) 
0.59 (140) 
0.64 (103) 


Full Band Fitting"* 
















Average 




347+g 


"■'t_o.o6 


456t^,^ 


7 n7+0.10 


230j:«« 


1.05 (327) 


Pole-on (Model 4 
Side-on 


in table ??)*= 


229t;f 
368ll?9 


f, 7,-1-0.07 
"•'J^-0.10 

"■'1-0.13 


34ltil? 


7gg+0.12 
'•""-0.15 
7 10+0.21 
'•^^■^-0.56 


282+164 
999+199 
^^^-105 


0.98 (243) 
0.82 (210) 



Pole-on phase; ip = 0.73 ± 0.25. Side-on phase; ip = 0.23 ± 0.25. 

The line center (keV). That of the fluorescent component is fixed at 6.40 keV. 

The determination of continumn spectrum is performed in 4.5 keV - 10.0 keV 

with a single A^h and single temparature model. 

The determination of continumn spectrum is performed in 0.8 keV - 10.0 keV 

with a multi A'^h and single temparature model (Model 4 in table ??). 

The temparature is fixed to the Ginga value. Continuum parameters arc fixed. 




Figure 17. The expected line beaming factor ^poio/^sido of V834 
Cen due to pure geometrical effects at the limit of large optical 
depth, calculated by considering the system geometry of V834 
Cen and the exposure for each phase. It is shown as a function of 
h/r. The solid line corresponds to the most likely system geometry 
of (i,/3) = (45°, 25°), while dashed lines reflect uncertainties in 
the system geometry. 



line, under the presence of the stronger He-like line adjacent 
to it. Furthermore, we expect the hydrogen-like iron line to 
be intrinsically less coUimated than the He-like resonance 
iron line, because the H-like line photons are produced pre- 
dominantly in the top regions of the accretion column (figure 
pi): there, the electron density is lower, thermal Doppler ef- 
fect is stronger, and the path of escape from the column is 
shorter, as compared to the bottom region where the He-like 
lines are mostly produced. Therefore, the H-like line is ex- 
pected to be less coUimated than the He-like resonance line 
(figure ns). Note in the observation blended He-like line is 
almost the same coUimation as H-like line. 

From these considerations, we conclude that the ASCA 
results on V834 Cen can be interpreted consistently by our 
POLE scenario. 



5.2 Determination of the accretion column 
parameters 

The line beaming effect we have discovered is expected to 
provide unique diagnostics of the accretion column of mag- 
netic WDs. The physical condition in the accretion column 
is described by four parameters: kT'^, h, r, and n°'\ To de- 
termine these parameters, four observational or theoretical 
constraints are required. Usually, observations provide two 
independent quantities, the temperature /cTqbs and the vol- 
ume emission measure VEM, which is ~ ('i-e'')^ • hr^. Also 
there is one theoretical constraint, that the shock heated 
plasma cools only by the free- free cooling, which relates 
fcT''', n^*", and h as in equation (A5). 

With these three constraints, and taking into account 
the Aizu's solution (1973), we can express h and r as 



h 



1.61 X lO^cm 



2.15 X 10 



7 f 
cm — 

VI 



17.7keV/ \^10«cm 
VEM 



(12) 



30 X 1054cm 



-3 



17.7keV J 



lOiScm" 



(13) 



Here, we normalized the VEM to the value of 1.30 x 10^* 
cm~^ obtained from V834 Cen, adopting a distance of 86 
pc (Warner 1987). The value of kT""" = 17.7 keV was deter- 
mined from the expected mass of WD by the observed ratio 
of H-like to He-like iron K^ lines, considering the vertical 
temperature gradient (Ezuka and Ishida 1999 figure 5); it 
is consistent with an independent calculation by Wu et al 
(1995) based on the Ginga observation (Ishida 1991). When 
nl'^ is low, the solutions to h and r imply a long cylinder- 
like column, while a flat coin-shaped geometry is indicated 
by high values of n"^. However, due to the lack of one more 
piece of information, the value of n^'' has so far been left 
undetermined. As a consequence, we have not been able to 
determine the column geometry. 



© 2001 RAS, MNRAS 000, Q-^ 



X-Ray Beaming due to Resonance Scattering in the Accretion Column of Magnetic Cataclysmic Variables 13 




1e+14 



1e+15 1e+ie 1e+17 

dtii^il.> [cm-''] 



Figure 18. The expected beaming factor ^poic/faide of V834 Cen, 
calculated by Monte-Carlo simulations under the actual observing 
condition with ASCA, and shown as a function of the assumed 
post-shock density. The shock temperature is set to 17.7 keV, and 



the volume emission measure is set to 1.30 X 10 cm 



The 



dashed and dotted data points indicate ^ of He-like and H-like 
resonance iron K^ photons, respectively. The solid data points 
show 5 of the blended He-like iron line, which includes intercom- 
bination and forbidden lines. The error bars take into account 
those of the geometrical parameters. 



The present work provides us with the needed fourth in- 
formation, the value of 5, which reflects the accretion column 
condition. Suppose we specify a value of nf'. Then, equations 
(h2) and ( |l3| ) determine h and r respectively, which in turn 
are used as inputs to the Monte-Carlo simulation to predict 
^. In this way, we have calculated ^ for V834 Cen as a func- 
tion of rfj^, and show the results in Figure hq. By comparing 
it with the observed value of (^poio/Caido )"'''' ~ 1.87, we ob- 
tain the best estimate as 72^'' ~ 4.6 x 10^* cm~"^. This in turn 
fixes the column height as ft ~ 3.5 x 10^ cm and the radius 
as r ~ 3.2 x 10^ cm. These values are considered typical for 
polars. 

Currently, the errors on ^poio/C=idc are so large that the 
value of n°'' is uncertain almost by an order of magnitude, 



with the allowed range being 1 x 10 



1 X 10' 



However, our new method will provide a powerful tool for 
next generation instruments with a larger effective area and 
an improved energy resolution. 



of metal abundances of polars, 0.1 - 0.8 solar, measured by 
Ezuka & Ishida (1999) is considered still valid. 

5.4 Mystery of POLEs 

Our Monte-Carlo simulation (§H) and the observed study of 
V834 Cen (SMJ) consistently indicate that the resonance iron 
K lines are enhanced by a factor of ^„ = 2 - 2.5 in the axial 
direction of accretion column. The effect is thought to be 
ubiquitous among polars, because we have so far employed 
very typical conditions among them. Consequently, if a polar 
with co-aligned magnetic axis (/3 ~ 0) is viewed from rear 
the pole-on direction (z ~ 0), we expect the iron K line 
EW to be persistently enhanced by 2 - 2.5 times. When 
this enhancement is not considered, such objects would yield 
artificially higher metallicity by similar factors. We conclude 
that the three POLEs (9l[) are exactly such objects. 

In condition oi ( > 2, the extremely high face- value 
abundances of AX J2315— 0592 will be modified to be ~ 1 so- 
lar, and abundance of RX J1802. 1-1-1804 or AX J1842-0423 
comes into the measured distribution of other polars within 
error. We can expect a strong coUimation with an adequate 
density of ria ~ 10^® cm~^ (figure |lO|b) and an adequate 
shape of h/r ~ 0.1-0.5 (figure hda) with typical radius ~ 10^ 
cm. Thus, one example of strong coUimation calculated for 
He-like iron resonance line is in the condition that temper- 
ature of lOkeV, r = 5 X lO'' cm, n, = 7.9 x lO^'^ cm"^ 
and VEM = 10^^ cm~'^, which yield an X-ray luminosity of 
1.0 X 10^^ erg s"^ (2 - 10 keV). Hereafter, we call this condi- 
tion the strong case. Figure Il9 shows the expected enhance- 
ment (" of iron He-like line (sum of resonance line, intercom- 
bination line and forbidden line) with various geometrical 
conditions (i and (3). In figure EOl we have converted the 
result of figure M into cumulative probability distribution 
of (". We hence expect (" > 2 at the condition of strong case 
shown above, with a probability of ~ 4.5 %. This estimate 
is in an rough agreement with the observation, i.e. the three 
POLEs among the known ~ 50 polars. We reconfirms that 
the iron abundances derived from the observation (Ezuka 
& Ishida 1999) remains valid to within 30 % for a major 
fraction of polars. Thus, we conclude that the mystery of 
POLEs have been solved by the proposed line-coUimation 
mechanism. 



5.3 Effects on the abundance estimates 

The iron abundance of other polars, measured by Ezuka & 
Ishida (1999), are subject to some changes when we properly 
consider the resonance scattering effects. For this purpose, 
we show in table H the geometrical parameters (i and /3) of 17 
polars, which have been randomly selected from currently- 
known ~ 50 polars. The expected enhancement C = X/ h • ^ 
are also listed with assumption that the resonance lines are 
collimated as much as that in the case of V834 Cen. Thus, 
the implied corrections to the iron abundances of these ob- 
jects are at most ±30 %, which is generally within the mea- 
surement error. Polars with almost side-on or almost pole-on 
geometry, like POLEs, are subject to a relatively large mod- 
ification in the abundance estimates, but such a geometry is 
a rare case (here, only VY For) . Therefore, the distribution 



6 CONCLUSION 

In order to explain the extremely intense iron K lines ob- 
served from several Galactic X-ray sources (including two 
polars; section hi), we have developed a scenario of POLEs, 
that the iron lines from accretion poles of magnetic WDs are 
axially collimated by the two mechanisms (section H) which 
become operational under large optical depths for the reso- 
nant line scattering. One is the geometrical effect which be- 
comes effective when the accretion column is rather short, 
while the other is a physical mechanism that the Doppler 
shifts due to vertical velocity gradient in the post-shock flow 
reduces the cross section for resonance scattering along the 
field lines. 

In section H, we have carried out Monte-Carlo simula- 
tions, and confirmed that the velocity-gradient effect, aug- 
mented by the geometrical effect, can enhance the iron K 
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Table 4. Geometrical parameters of randomly sampled polars, and values of C expected for them 



object name 


i" 


/3" 


c 


reference 


BLHyi 


70 ±10 


153 ± 10 


0.70 


Cropper 1990 


UZ For 


-88 


~ 14 


0.71 


Ferrario et al. 1989 


VVPup 


76 ±6 


152 ±6 


0.71 


Cropper 1990 


ST LMi 


64 ±5 


141 ±4 


0.71 


Cropper 1990 


AN UMa 


65 ±20 


20 ±5 


0.87 


Cropper 1990 


QQ Vul 


65 ±7 


~23 


0.88 


Cropper 1990, Schwope et al. 2000, Catalan et al. 1999 


V1309 Ori 


~80 


~30 


0.86 


Harrop-AUin et al. 1997 


EFEri 


58 ±12 


27 ±18 


0.96 


Cropper 1990 


ARUMa 


50 ±10 


23 ±13 


1.05 


Szkody et al. 1999 


WWHor 


-74 


-48 


1.06 


Bailey et al. 1988 


DP Leo 


76 ±10 


103 ±5 


1.06 


Cropper 1990 


MRSor 


43 ±5 


38 ±5 


1.09 


Cropper 1990 


J1015+0904 


55 ±5 


43 ±7 


1.10 


Burwitz et al. 1998 


V834 Cen 


45 ±9 


25 ±5 


1.11 


Cropper 1990 


AM Her 


52 ±5 


66 ±5 


1.13 


Wickramasinghe et al. 1991, Isliida et al. 1997 


EK UMa 


56 ±19 


56 ±19 


1.14 


Cropper 1990 


VYFor 


9±3 


-9 


1.66 


Beuermann et al. 1989 



Inclination angle (degrees). 

Pole colatitude (degrees). 

Expected beaming factor to the average flux, assuming 

that the same collimation as the case of V834 Cen occurs 

and that only single accretion column omitts. 




10 20 



40 50 



Figure 19. The expected enhancement (^ to the average line intensity calculated on the plane of {i,l3), assuming that accretion column 
exists only on one pole with the same collimation as the case of V834 Cen (left) or the strong case shown in text (right). Contour levels 
are shown in the figure. The polars listed in table hi are also plotted in the left panel with tick marks. 



line EW in the pole-on direction up to a factor of 3.0 as 
compared to the angular average. This is higher than the 
maximum collimation available with the geometrical efTect 
alone (Appendix N) . 

In order to experimentally verify our interpretation, we 
have analyzed in section the ASCA data of V834 Cen, 
which has a suitable geometry in that our line-of-sight to 
the accretion column changes from 20° to 70° as the WD ro- 
tates. Through detailed phase-resolved X-ray spectroscopy, 
the EW of the He-like iron-Ka hne has been confirmed to be 
enhanced by a factor of (^poio/^sido)"'"' = 1.87 ± 0.54. In sec- 
tion 5.1, we have examined whether the observation can be 



explained away by any other mechanism, and we have con- 



cluded that this observational result of V834 Cen strongly 
reinforces our interpretation of POLEs. 

Although the resonance lines are colhmated with the 
proposed beaming effect, the previous measurements of the 
distribution of metal abundances of polars is considered still 
valid (section 5.3) except for POLEs. With proposed mech- 
anism, the extremely high face- value abundances observed 
in POLEs can be reconciled with the average abundance 
measured from the other polars. Thus, the POLE scenario 
successfully solves the mystery of the extremely strong iron 
lines observed from the three X-ray sources. 

In addition, our scenario provides a new method of 
unique determination of physical condition in the accretion 
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Figure 20. The cumulative probability distribution of finding 
f higher than the specified threshold fth, assuming that objects 
are randomly distributed in terms of i and (3. The dotted and 
solid curves represent the case of V834 Cen and the strong case, 
respectively. 



column, using the beaming factor of r eson ance lines as a 
new observational information (section 5.2). This will be a 
powerful method for the next generation instruments. 



We thank the members of the ASCA team for spacecraft 
operation and data acquisition. 



APPENDIX A: PLASMA IN AN ACCRETION 
COLUMN OF A WD 

A flow of accreting matter captured by the magnetic field 
of WD channels into an accretion column, where a standing 
shock forms and heats up the matter to a temperature kT"^ 
of 

3 GMv/nfiniH 

8 



kT" 



= 16 



(o 






615 



O.5M0 



(R^\ keV, (Al) 
V 109cm/ ' ^ ^ 



where k is the Boltzmann constant, G is the gravitational 
constant, AfwD is the WD mass (typically 0.5Mq), and 7?wd 
is the radius of the WD (typically 10^ cm). So the plasma 
has a typical temperature of hard X-ray emitter, forming 
an accretion column illustrated in figure hi The velocity be- 
neath the shock front u"^ is described with relation to the 
free-fall velocity itff as 



Mff 

4 



= 0.9 X 10" 



2G'MwE 

/tWD 

MwD_ 





0.5M( 



1/2 



r?v 



V lO^cm 



-1/2 



i/s. (A2) 



Assuming that the plasma is a single fluid and that the 
abundance is one solar, the electron density of the post- 
shock plasma is given as 



M 



nr-^u"^fj,mu 



7.7 X 10' 



X 0.518 

M 



IQlBg s-1 



A/v 



O.5A/0 



-1/2 

(.ftWD \ 
lO^cm/ 



1/2 



VS X lO^cmJ 



(A3) 



where the value 0.518 is the fraction of electron density to 
the total density assuming the solar abundances. In the ac- 
cretion column, kT, u and Ue. all have a vertical gradient 
from the shock front toward the WD surface. Numerically, 
the vertical profiles of these quantities as a function of the 
distance z from the WD surface, normalized by h, are cal- 
culated by Aizu (1973) as 



kT 
kT"" 






(A4) 



where each quantity is normalized to its value immediately 
below the shock front; fcT"'', u"'^ , and n"^ . 

The column radius of the accretion column is typically 
r = 5 X 10^ cm. Since the shock front is sustained by the 
pressure of the post-shock plasma against the gravity, h is 
described by free-free cooling time scale itf.cooi of the heated 
plasma as ft ~ u°''tfi_oooi- According to Aizu (1973), h is given 
more specifically as 



h 



0.605?i"''iff.oooi 

^7 f kT'"^ 



= 1.9 X 10' 



\^ 16keV 



7.7 X 1015cm" 



2m(A5) 



where ttt.cooi is given by 
3nefcT='' 



Ut., 



0.35 



26ff 



ll6keVy 1,7.7 x1015 cm 



s (A6) 



with Eft being the volume emissivity of free-free emission [eq. 
(5.15) in Rybicki & Lightman 1979]. Note that the ion to 
electron energy transfer time scale tcq 






5.6 X 10" 



16keV 



7.7 X 1015cm 



[see eq (5.31) in Spitzer 1962] is much shorter than ttt.cooi, 
so the ions and electrons are thought to share the same 
temperature. 



At the temperature of a few tens keV [equation (Al)], 
the electron scattering dominates the opacity in the hard X- 
ray band. Actually, the optical depth of the accretion column 
for free-free absorption, Ta, is given, relative to the electron 
scattering optical depth tt [equation (|l|)] as 



T{{ 



0.88 X 10" 



Tt 



(-^y (■ 



kT 



v 



V6keVy V16kcVy 



7.7 X 1015cm 



(A8) 



[see equation (5.18) in Rybicki & Lightman 1979]. Thus, 
the free-free absorption is negligible compared to Thomson 
scattering. 

The cross section of resonance scattering ctrs for a pho- 
ton with energy Eo can be described generally as 



Tre i 

crs = ./i2 , — ■ — 

m-iC V27rA£; 



exp 



{Eo — Ejii 
2A£2 



(A9) 



where /12 is the oscillator strength for Lyman-a transition 
(n = 1 to 2), Eb.s is the resonance energy in the rest frame. 
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Figure Bl. The coUimation factor of geome tric al beaming. The 
flux at S = to the averaged flux [equation (p2|)] . 



and /S.E is a resonance energy width [eq. (10.70) in Rybicki 
& Lightman 1979]. Numerically, the first factor is -^ — /12 = 
1.2 X 10^^^ cm~'^, and the energy width /S.E in the second 
factor of Gaussian is determined by equations (g) and (|6|). 
Then, the cross section of resonance scattering at the line- 
center energy is given as 

ans = 2.0 X 10-^« (^7^1^) f T7^) '"'^ cm^ (AlO) 



6.695kcVy Viekevy 



APPENDIX B: GEOMETRICAL BEAMING IN 
THE ACCRETION COLUMN 

How much enhancement can we expect in "geometrical 
beaming" in the accretion column (see §2.1)? Consider the 
case that the resonance line photons can only escape from 
the surface of the column. The directional photon flux 
emerging from the column is given as 



f{.d)d 



cos 9 (X nr cos 6 + 2rh sin ( 



(Bl) 



where 6 is the angle measured from the column axis. There- 
fore, the flux along S ~ is enhanced by a factor 






r2-f- fr/l 



1 + 1^ 



(B2) 



where < / > is the average of f{9) over 9. At the coin-shaped 
limit (-^ -^ 0), this factor approaches 2 (figure Bl). 
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